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Abstract. 

A new singular perturbation method based on the Lie symmetry group is presented 
to a system of difference equations. This method yields consistent derivation of a 
renormalization group equation which gives an asymptotic solution of the difference 
equation. The renormalization group equation is a Lie differential equation of a Lie 
group which leaves the system approximately invariant. For a 2-D symplectic map, 
the renormalization group equation becomes a Hamiltonian system and a long-time 
behaviour of the symplectic map is described by the Hamiltonian. We study the 
Poincare-Birkoff bifurcation in the 2-D symplectic map by means of the Hamiltonian 
and give a condition for the bifurcation. 

PACS numbers: 02.20-a, 02.30Mv, 02.30Oz, 02.40Xx, 45.fOHj 
1. Introduction 

There have been many studies concerning application of renormalization group method 
of quantum field theory as a singular perturbation method to treat differential equations 
since the work of Illinois groups [T][2]- Although this method enables us to remove 
secular or divergent terms appearing in a naive perturbation solution by renormalizing 
integral constants appearing in the lowest-order of perturbation solution and to give 
a well-behaved asymptotic solution, it is necessary to calculate a naive perturbation 
solution to derive the renormalization group equation, and whether the system could 
be renormalized or not depends on the functional form of a naive perturbation solution. 
The present authors presented another renormalization group method in terms of 
the Lie symmetry group [3]. Without calculating a naive perturbation solution, the 
renormalized solution is constructed as an invariant solution under transformation which 
leaves the system approximately invariant. The approximate Lie group transforms an 
unperturbed solution into a renormalized perturbed solution. The both renormalization 
group methods were based on the theory of continuous group and so they were developed 
to apply a system of differential equations. For some discrete dynamical systems such 
as a system of difference equations, there has been found a working procedure to derive 
a renormalized difference equation by means of a naive perturbation solution [1][5]. 
However, the concept of group was not used in the procedure to derive a renormalized 
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discrete system and so legitimacy of the result is not clear. The purpose of this paper 
is to develop a renormalization group method in terms of the Lie symmetry group 
for a system of difference equations. Apparently, it seems to be difficult to use the 
same procedure as in the case of a differential equation when we treat a discrete system 
because Lie symmetry groups are one-parameter continuous groups. However, we regard 
a discrete system as an algebraic system and consider a particular class of approximate 
Lie symmetry groups which transforms a perturbed system to an unperturbed system 
approximately. Then we succeed in obtaining a renormalized perturbed solution by 
means of the so-called Lie differential equation corresponding to the approximate Lie 
symmetry. Thus, we present a renormalization group method for a system of difference 
equations, where the continuous group theory is fully employed, and the Lie equation 
becomes the renormalization group equation for a difference system. As an application 
of the new method, we study the Poincare-Birkoff bifurcation in a 2-D symplectic map 
and give a condition for the bifurcation explicitly. 



2. Renormalization group method with Lie symmetry 

Let us consider the following 2-D symplectic map of action-angle type, (u n , v n ) i— > 
(u n+ i,v n+ i): 

u n+l = u n + V n+ \, 

v n+1 = v n + au n + eg(u n ), (1) 

where the coefficient a G M is constant, g is a function of u n , and e is a perturbation 
parameter which is small. The system (CQ) reads 

(2) 



(3) 
(4) 





Under a transformation of coordinates (u n ,v n 

• exp(iu;/2) 




2cos(u;/2) 2sino)coso; 

1 ■ cxp (-icv/2) 

2cos(u;/2) 2sinojcosoJ 





(5) 



exp (— itu/2) exp(ia;/2) 
-2isin(w/2) 2isin(u;/2) 

the linear part of ([2]) is diagonalized as follows: 

z n +i \( e-' luJ \ f z n \ ^ g(u n ) ( iexp(-io;/2) 
z n+ i J I e luJ I \ I 2 sin u I — iexp(zu;/2) 

where a + 2 =: 2cosu;, and denotes the complex conjugate of z n . Here < a < — 4 
because we are interested in a case the origin is assumed to be elliptic. 

Since the second component of ([5]) is the complex conjugate of the first, we have 
only to consider the first component of the equation, which is 

:, , iexp(— ia;/2) . , . . 

z n+1 = e^zn + e ^ —g{u n ). 6 

2sinw 
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Let admit a Lie group transformation whose infinitesimal generator takes the 

form 

X(n,z n ,z^,z n+ i,z^i) = d £ + r] z (n,z n ,z^)d Zn + rf{n, z n , X)c\^ 

+r] z (n + 1, z n+1 ,z^)d Zn+1 + r] z (n + 1, z n+1 , z n ^)d^ TT . (7) 

We can rewrite (H) to 



X(n, z n , Zn, z n+1 , z n+1 ) =d £ + r]{n, z n , z n )d Zn + r](n, z n , z n )c\^ 



+r](n + 1, z n+1 , z n+1 )d Zn+1 + r)(n + 1, z n+1 , z^d^, (8) 

because it can be shown that rf = rf. Then the determining equation for ([6]), which 
determines r](n,z n ,~z^) of the vector field (IE]) |6J, is given by 



. iexp (— iuj/2) s 

X [n, z n , z n , z n+ \, z n+ i) < z n+ i — e z n — e — 9\U n/ 

2 sin uj 



= 0. (9) 

Equation \Q\ 



Because we wish to find such a symmetry that leaves the system approximately invariant 
to leading order, we need only to solve the following leading-order determining equation: 

. . f iexp(— \uj/2) y 

X (n, z n , z n , z n+ \, z n+ i) < z n+ i — e z n — e — 9\u n/ 

2smu 



= 0(e), (10) 

z n j r \ — e z n 



r](n,z n (e),z n (e)). (13) 



^=>r](n+l,e w z n ,e^z n ) -e ""77(71, z n , z n ) = 16X ^, ™^ g{u n ). (11) 

2 sincj 

By solving (TIT]) , we obtain the infinitesimal generator X admitted by the system ([6]) in 
the leading-order approximation. 

Using the obtained infinitesimal generator X, we construct a group-invariant 
solution of the system ([6]). The group- invariant solution, z n = z n (e), satisfies the 
following Lie equation: 

X{zn-z n (e)}\ Zn=Zn{£) =0, (12) 

which reads 

dz n (e) 
de 

We refer to the Lie equation (1131) as the renormalization group equation. Solving 
the renormalization group equation by adopting solutions of unperturbed system as 
a boundary condition, i.e. 

z n (e = 0) = z<®, (14) 

where z n = z^ denotes the solution of the unperturbed system, we obtain an asymptotic 
solution of the system. 

3. Power series nonlinear terms 

Let us consider a case g{u n ) is a power series with respect to u n , i.e. we set 

oo 

g( u n) :=^A jU i. (15) 

3=0 
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where Aj G R. Under the linear transformation (jJJ), (1151) can be written as 

g(u n (z n ,z^)) = S ^Bi„ l (uj)z l n z^ n , (16) 

l,m 

where 

B im {u) = A l+m i +m Ci exp [-i(Z - m)w/2]. (17) 
Here := fe!/[r!(A; — r)!]. Then the determining equation (TTTj) becomes 

^(ra + 1, e~ lw z n , e 1<J ^) - e~ lu) n(n, z n ,z^) = iex P(^ W 2 ) ^(q;)^^". (18) 

Z Sill w 

The nature of the solution of ([TBI depends on whether u;/27r is rational or irrational. 
3.1. uj/2tt is irrational 

Let uj/2ti be irrational, then only terms proportional to z^f 1 ^ 1 cause resonance, and 
the solution of the determining equation ([TBI becomes 

iexp i\uj/2)Bi m (uj) 



r](n,z n ,z n ) = ^2 



2 sinu; 

l,m s.t. !— m=l 



+ y- iexp(-iu}/2)B lm (u) 

/ -— ' 2 sin u; iexp I— i(Z — m)u/| — exp [— icul } n n 

l,m s.t. i-m^l 7 J * 1 J J 

Here exp [— i(Z — m)a;] — exp [— ia;] 7^ for I — m 7^ 1. The terms proportional to n are 
called secular terms. The renormalization group equation i.e. ( TT3l) becomes 

d2„ _ ^ iexp (iuj/2)Bi m (uj) l 



de *—f 2 sin u; 

t,m s.t. l—m=l 



+ iexp (h;/2)B lm (u) 

' 2 sinu; jexp [— i(l — m)uj] — exp [— iu]} n n 

l,m s.t. l-m^l ^ ' J ^ 1 J J 

Because we are interested in a long-time behaviour of the system, we consider the case 
of n ^> 1. Then we neglect non-secular terms and the renormalization group equation 
becomes 



dz n iexp(iu>/2) ^ 2l 

fif / sin m < * 



de 2 sin uj 

1=0 



„ / A 2 l + l 2l+lCl n \ z n\ 21 Z n , (21) 

2 sin u; z — ' 

where we use ( ITTjl . Because d|z n |/de = 0, introducing R := |,z n |, the solution of (T2TT) 
becomes 

= z n (0) exp [ — — V A 2l+1 2l+1 dR 2l ne ) . (22) 
\ zsinw *— ' / 

Thus, we obtain an asymptotic behaviour of a difference equation by solving a differential 
equation. 
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3.2. lo/2h is rational 

Let a;/27r be rational i.e. u = 2 r Kqjp where q,p e N, then denominators of the second 
terms in the right hand side of ( 120|) . exp [— i(Z — m)uj] — exp [— iu], become zero when 
I — m = 1 + kp for k 6 Z. Therefore terms z l n z n m where I — m = 1 + kp give secular 
contribution to the determining equation ( fl8l) . The solution of the determining equation 
becomes 

iexp {uo/2)Bi m (uj) 



r)(n, z n , z n ) = ^ 

l,m s.t. l—m=l+kp 



~^Z n Z n 



+ 



E 



2 sintu 

iexp (-io;/2)S Zm (a;) 



2 sin to {exp [— i(Z — m)a>] — exp [—ico] } 



l,m s.t. I— m^l+kp 

In the case of n ^> 1, the renormalization group equation reads 
dz n x v iexp (iu;/2)5 /m (^) ; _ m 



de 



E 



Z,m s.t. (— m=l+fcp 

which is rewritten as 

dz n iexp (ko/2)n 



2 sin a; 



(23) 



(24) 



d£ 



2 sinu; 



E* 

Z>l,m>0 



Z+mp Z-H^Jp? 



|2(Z~1) mp+1 



i+mp- 1 ('- L I | Zn ^ 



Z>0,m>l 



(25) 



or 



dz n iexp(ia;/2) 
d7 = 



2 sin a; 
i 



./>l,m>0 



r 



|2(Z-1) mp+ 



E« 



Z Z+mp-H^J |^n| & 



21^-mp— 1 



Z>0,m>l 



2 sincj 



^] ^-21+mp — 1 2Z+mp— 1 +mp | I 



2(/-l) r mp+l 
n 



|2Z— -mp-1 



4>l,m>0 

+ ^ ] ^21+mp-l 2l+mp-lCl\z n \ M ' Z n 
Z>0,m>l 

where r := en. This is a Hamiltonian system whose Hamiltonian is 



H(z nj z Tl 



E 



l > 0,m > 
(Z, m) ^ (0, 0) 



L4 



2Z+mp-l" 



(2l + mp-l)\ 
(l + mp)\l\ 



21 ( z mp _|_ —mp\ 



(26) 



(27) 



Because (1261) describes also n dependence of z n (e), the phase space structure of this 
Hamiltonian system provides a long-time behaviour of the difference equation. 



3.3. u}/2tc is close to a rational number 

In order to study the Poincare-Birkoff bifurcation in the symplectic map, let us consider 
the case lv/27t is close to a rational number, 

uj = uj q + e5, (28) 
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where uj /2it is a rational number, ujq = 2irq/p. The difference equation (jH]) becomes 



Zn+l 



-UJQ 



z n + e 



iexp (— ia>o/2) 



-g[u r 



+ 0(e 2 ). 



(29) 



2 sin uq 

With the same procedure as we have followed in the previous subsections, we can 
derive the renormalization group equation in the case of n ^> 1 as follows: 



dZr, 



dr 2 sin u 



£ -4 



21+mp—l 2l+mp—lClJ rm p\Z n 



\2{l-l) mp+X 

Am 



4>l,m>0 

I \^ A ry I _ 12/^- 

~r / sT-21+mp-l 2l+mp-l i ~'l\ z n\ z n 

i>0,m>l 

For simplicity, let us consider such a case as follows: 
:= 

p-i 



2l——mp—l 



iSz n . (30) 



6 J2 p-i C i ex P ["i(2/ - V + l)^/2]z l n l 



-i-i 



(31) 



1=0 



where b is a parameter. Then the renormalization group equation (|30|) , which describes 
a long-time behaviour, becomes 



dz„(T) 



ib 



OL- 



i6 



z^f i5z n , 



(32) 



dr 2sinu;o P ± 2 2sina;o 

where a = 1 when p is an even and a = when p is an odd. This is a Hamiltonian 
system whose Hamiltonian is 

ib ,„ ife 

n_l C _ 

2 



H = a- 



p sin u; 
At fixed points z n satisfies 
ib 

2 sin u;q 



a- 



+ 



2p sin a; 
ib 



(33) 



"3r 



■p-1 



Setting 2 n := re 16 * where r and 6* G M 
sin(p#) = 0, and r p ~ 2 



2 sin ujq 

25 sin oft 



i5z„ = 0. 



> 0, 



b ( a p __iCV + cos(p6) 



which reads 



25 sin uiq 



„p-2 = 

6 ^a p _iCE + 1 
„ 2 25 sin ujq 



for 9 



TV 



21. 



P 



71 



b ( a v ^\Cv — 1 



for 6= -(21 + 1), 
p 



(34) 
(35) 

(36) 
(37) 



where I = 0, 1, ... ,p — 1. When p is even number which satisfies p > 4 and 65 > 0, 
there are 2p fixed points. In this case, it can be shown that the half of them are elliptic, 
while the others are hyperbolic, and the Poincare-Birkoff bifurcation occurs. Here the 
resonance structure in this case consists of a chain of p resonant islands. 
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4. Concluding remarks 

The renormalization group method with Lie symmetry is developed and applied 
to difference equations such as symplectic map. When we apply the conventional 
renormalization group method to a singular perturbation problem of difference 
equations, it is difficult to find how to construct the renormalization group equation, 
which is a continuous differential equation. In this paper, we succeed in the consistent 
derivation of a renormalization group equation to a difference equation by employing the 
theory of Lie symmetry. Our results are significant in the sense that a long-time behavior 
of a perturbed difference equation is described by a continuous differential equation, that 
is the renormalization group equation. Furthermore, for the case of a simplistic map, it 
is shown the renormalization group equation becomes a Hamiltonian system. Then the 
phase space structure can be described by the Hamiltonian. As an application of our 
theory, we analyze the Poincare-Birkoff bifurcation in a 2-D symplectic map and the 
condition of the Poincare-Birkoff bifurcation is given to the 2-D symplectic map. 
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